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ABSTRACT 

The late-time phase-ordering kinetics of the 0{n) model for a non-conserved order 
parameter are considered for the case where the 0{n) symmetry is broken by the 
initial conditions or by an external field. An approximate theoretical approach, 
based on a 'gaussian closure' scheme, is developed, and results are obtained for the 
time-dependence of the mean order parameter, the pair correlation function, the 
autocorrelation function, and the density of topological defects [e.g. domain walls 
{n = 1), or vortices {n = 2)]. The results are in qualitative agreement with exper- 
iments on nematic films and related numerical simulations on the two-dimensional 
XY model with biased initial conditions. 
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1 INTRODUCTION 

The field of pliase-ordering kinetics lias seen a resurgence of interest in recent years 
It is by now well established that in the late stages of ordering a scaling regime 
is entered, characterised by a single time-dependent length scale L{t). The exis- 
tence of a single characteristic scale has the consequence, for example, that the 
pair correlation function for the order parameter field exhibits the scaling form 
C{r,t) = f[r/L{t)], where f{x) is a scaling function. Systems with nonscalar order 
parameters have attracted special attention recently, partly in response to exper- 
imental interest in liquid crystal systems p|. In particular, approximate methods 
developed for calculating scaling functions for scalar fields P, || , have been extended 
to vector fields |], |^. So far, however, little effort has been devoted to the role of 
external fields, or a symmetry-breaking bias in the initial conditions. 

In this paper we study the phase ordering dynamics of systems with non-conserved 
order parameter, in situations where an external field is present and/or the initial 
state contains a bias. We may call these systems 'off-critical', as opposed to critical 
systems whose initial state and dynamics are rotationally invariant (or inversion 
symmetric for scalar fields). An immediate difference relative to critical quenches is 
that the system will (in most cases) reach global saturation exponentially fast, thus 
limiting the duration of the scaling regime. We study both scalar and vector sys- 
tems. In the former there is a simple magnitude relaxation of the order parameter, 
while in the latter there is also a slower orientational process. Physical realizations 
are spin systems (magnets), nematic liquid crystals, and order-disorder transitions 
in binary alloys. The problem is addressed using an analytically solvable approach, 
the systematic approach (SA) which is a form of 'gaussian closure' approxima- 
tion. At leading order, it reproduces the results of earlier approximation schemes 
0, ^, but has the virtue that it is (at least in principle) systematically improv- 
able [^. We derive an approximate linear equation for an auxiliary field m. This 
equation is exact in the limits — * oo or d oo, where n and d are the spin and 
space dimensions respectively. From this we obtain predictions for the mean order 
parameter, the topological defect density, and the pair correlation function. 

The systematic approach was developed by Bray and Humayun (BH) [§] for 
simple non-conserved dynamics with no field or bias, in which case the results of 
Ohta, Jasnow and Kawasaki (OJK) [Q], and the extension thereof to vector fields 
by Bray, Puri and Toyoki (BPT) |^ are recovered in leading order. This approach 
has the virtue, however, that it can in principle be systematically improved. We are 
not concerned here with the possibility of improving the leading order calculations, 
which is a problem of great technical difficulty |]TD[ , but in showing that the approach 
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(in lowest order) can be successfully extended to off-critical systems. 

In section ^ we discuss the original systematic approach. For simplicity we begin 
with scalar fields. The extension to vector fields is straightforward. In section ^ we 
generalize the scheme to include a time-dependent external field and an initial bias. 
In this case the relevant phenomenon is different for scalar and vector systems and 
so the construction of the approach is also different. For scalars the important effect 
of the field is a domain wall driving force, while for vectors the bulk rotation of 
the order parameter plays an important role. The magnetization is compared with 
an earlier prediction from the large-n solution of the time-dependent Ginzburg- 
Landau (TDGL) equation, by Bray and Kissner (BK) We also compare the 

predictions for the planar XY model with an initial bias with data from experiment 
and simulations on nematic liquid crystals [|l^ . 



2 Critical Quenches 

An appropriate free-energy functional to describe a phase ordering system is the 
Ginzburg-Landau form 



F[$]= Jd'x (V0)' + V(0) 



(1) 



where V{(f)) is a Mexican hat potential, e.g. V{(j)) = (1 — 0^)^, with a degenerate 
ground state manifold 0^ = 1. For scalar systems, the potential has the usual 
double-well structure, and the ground state has a discrete, two-fold degeneracy. 

An appropriate equation of motion for non-conserved dynamics is given by the 
TDGL equation 

d(p SF .7 dV 

TT- = ^ = ^ . (2) 

We will consider scalar and vector systems in turn. 

2.1 Scalar fields 



According to the seminal work of Allen and Cahn [|13|, the motion of the interfaces 
in non-conserved systems is determined, in the absence of a bulk driving force, solely 
by their local curvature. The detailed form of the potential V{(f)) and the particu- 
lar distribution of initial conditions are not important to the late-stage dynamics, 
and therefore to the scaling properties. The only relevant feature of V{(j)) is the 
double-well structure, which generates and maintains well-defined interfaces, whilst 
its detailed shape only affects the short- distance behaviour, such as the domain-wall 
profile. Similarly, the details of the initial random configuration determine the early- 
stage locations of the walls which, in the absence of any systematic bias, should be 
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irrelevant to the scaling properties. The systematic approach exploits this freedom 
through a convenient choice of V{(j)) and initial conditions. 

The TDGL equation (Q) for the evolution of a non-conserved, scalar field (/)(r, t) 
reads (where primes indicate derivatives) 

^ = v^0-n0). (3) 

In the 'gaussian closure' schemes 0, 0, |14[ |, ^ a new field m(r, t) is introduced, 
which varies smoothly on the domain scale and whose zeros define the positions 



of the walls. Following Mazenko |14|, the transformation 0(m) is defined by the 



equilibrium planar profile of an isolated wall, which satisfies 

<P"{m) = V'{<P) , (4) 

with boundary conditions 

0(±cx)) = ±1 , 0(0) = . (5) 

With this choice, near a wall (where this can be regarded as fiat) m plays the role 
of a distance from the wall, whilst inside the domains m is typically of order L{t). 
Rewriting equation (|^) in terms of m, and using (H) to eliminate V, gives 

m = V^m - f 1 - (Vm?) . (6) 

(p'{m) ^ ' 

The non-linear profile function 0(m) depends on the potential V"(0), so for gen- 
eral potentials equation (^ offers no obvious simplification over the original TDGL 
equation (H). As argued above, however, the scaling functions should be insensitive 
to the potential, and thus to the details of the profile function. The results of the 
OJK 0] and Mazenko |14[ approaches, for instance, make no explicit reference to the 



potential. The key step in the present approach is to exploit this idea, by choosing 
a particular form of V^(0) such that equation (H) is simplified. 

Bray and Humayun (BH) chose the profile function to satisfy 

0"(m) = —m(f)'{m) , (7) 

which is equivalent, via (^), to a particular form of the potential, as shown below. 
Integrating (|^ with the boundary conditions (|^) gives the wall profile function 

(f){m) = ^/2/7f dx exp(^-x'^/2^ 

= erf(m/V2) , (8) 
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where erf(x) is the error function. On the other hand, integrating (|^) once, with 
the zero of the potential defined by V{±1) = 0, gives 

V{<P) = (1/2) = (I/tt) exp(-m2) (9) 

= (l/7r)exp(-2{erf-i(0)}') , (10) 

where erf~^(a;) is the inverse error function. Note that @ is a monotonic mapping, 
with the constraint 0^ < 1 due to the boundary conditions (^. It follows that (|TOD 
only determines V{(j)) within the region 0^ < 1. For instances, we have V{(f)) ~ 
l/vr - 02/2 for 0^ ^ ^hile 1/(0) ~ (1/4) (1 - 0^)^ | ln(l - 02)| for 0^ ~ 1. At 
a temperature of T = 0, however, if the initial condition satisfies 02(r,O) < 1 
everywhere, the equation of motion (^ implies that (f)'^{r,t) < 1 everywhere at any 
later time. So 0(r, t) does not depend on V{(j)) for 0^ > 1, and there is no need to 
know the potential in this region (if T > 0, however, the points = ±1 must be 
global minima of V{(j)), in order to ensure stability against thermal fluctuations). 
With the choice (0), equation (P) reduces to 



m 



V^m + (l - (Vm)2) m . (11) 



This equation, though still non-linear, is much simpler than the original TDGL 
equation, while retaining the physical ingredients essential to describe the universal 
scahng properties. To make further progress, and to recover the usual OJK results, 
one simply replaces (Vm)^ by its average (over the initial conditions) 

{Vmf {{Vmf) . (12) 

Then ( pH]) becomes the self-consistent linear equation for m{r,t): 

m = V^m + a(t) m , (13) 

with 

a{t) = 1 - ((Vm)2) . (14) 

As in OJK's theory, to proceed with the calculations one takes the initial con- 
ditions for m to be gaussian distributed, with zero mean and correlator (in Fourier 
space) 

(mk(O)mk' (0)) = A{2TTY5{k + k') , (15) 



representing short-range spatial correlations at t = 0. Hence, from ([T3|) , m(r,t) is a 



gaussian field at all times. The formal solution of (0) in Fourier space is 

mk(t) = mk(0) exp (-kH + b{t)) , (16) 
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with the definition 

b{t) = f dt'ait') . (17) 







The function h{t) is determined self-consistently, as follows. Inserting (|l^)-(0) into 
definition (0) gives 

a{t) = db/dt = 1 - A exp(26(t)) J j^k'^ exp {^-2kH) . (18) 

Although this equation can be readily solved for the whole time range after the 
quench, it suffices for our purposes to extract the large-t behaviour of h{t) which is 
most easily obtained directly from (0). In the scaling limit we expect that a{t) 
(in fact we expect that a{t) ~ 1/t if all terms in (|T^ scale the same way). Neglecting 
the left-hand side of (|T^ and performing the momentum integral, gives 



with the definition 



d+2 

exp(6(t)) ~ (t/to)~ , t>to, (19) 



t^^ = — . (20) 

4(871)^^/2 ^ ^ 



From ([19|) we immediately get 



, , d+2 , , 

a{t) ^ , t > to , (21) 

confirming the consistency of the previous assumption. The short-time cut-off to 
was introduced as a device to prevent the breakdown of the momentum integral 
in ([18| ) as t ^ 0. It is equivalent to imposing a large-fc cut-off (of order l/tg ) 
associated with the smallest length in the system (e.g. the defect core size (), which 
removes the mathematical singularities. So to should be small. The condition t ^ tg 
characterizes the times for which the above asymptotic forms are valid, i.e. it defines 
the scaling regime {L{t) ^ () in which we are interested. Throughout this paper 
we will, however, for computational convenience, take t down to to while still using 
(p!9|). This implies that the calculated scaling functions will neither be accurate nor 
universal in the early-time regime t ~ to, but will not affect their behaviour in the 
true scaling regime. In this regime, the scaling functions will not depend on the 
specific way in which the cut-off has been introduced. 

Having determined 6(t), we may use ( pj] ) to write the explicit result for mk(t), 
valid for large-t, 

d+2 „ 

mk(t) = mk(0) (t/to) ^ exp{-kH) , (22) 
from which the two-time correlator, in Fourier and real space, follows immediately: 

(mk(ti)m_k(t2)) = A (tit2/ty ' exp[-e{U+t2)) , (23) 
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and 



^Vhh ( 4tit2 Y'^ ( r 



.2 



C„(1,2).M1M2)) = ^^^^J exp^-j^J. (24) 

where '1' and '2' are the usual shorthand for space-time points {^\^t\) and (r2,t2), 
and r = |ri — r2|. From this, we find the expected scahng form 

S,{t) ^ (m^) = ^ ~ L{tf , (25) 

in agreement with the physical interpretation of m as a length. It is interesting 
that the linear term a(t)m in (|T3|) , with a(t) given by (|2l|) , exactly reproduces the 
Oono-Puri extension of the OJK theory [|l^], designed to eliminate the implicit time- 
dependence of the interface-width, which is unphysical. The resulting difference in 
the time factors of (0) and (^), however, does not affect the scaling functions, 
which as usual only depend on the normalized correlator 

M1M2)) ^ ^ 

^(m(l)2) (m(2)2) V(*l+«2)V V 



Finally, we evaluate the correlation function of the field 0. Since from ( pSf ) m 
is typically of order at late-times, it follows from (|^) that the field saturates 
(0 = ±1) almost everywhere at late-times. As a consequence, we can make the 
usual simplification 0(m) = sgn(m) as far as scaling properties are concerned. The 
calculation of the average C(l, 2) = (sgn(m(l))sgn(m(2))) for a gaussian field m(r, t) 
is the same as in the OJK theory (see the Appendix). Therefore, the OJK scaling 
function C(12) = (2/7r) sin~^[7(l, 2)] is recovered. 

We conclude this subsection with a technical remark on the above procedure. 
To make the replacement (|12|) in a controlled way, BH systematize the treatment 
by attaching to the field m an internal 'color' index a, running from 1 to A^, and 
generalizing (|TID to 

rria = V'm, + - {XjN) J] (Vm^)' j m„ . (27) 

Equation (|ri|) is the case A^ = 1. In the limit N —>■ oo the 'color' sum is replaced 
by its average, giving the linear equation (0) (with m standing for one of the nia)- 
This procedure allows, in principle, for a systematic treatment of the results in 
powers of 1/A^. The replacement (0) is also justified in the limit d oo, when 
(Vm)^ is also a sum of a large number of random variables with relative fiuctuations 
~ \l\fd. This scheme, though less simple to systematize 0], makes clear that the 
leading order results become exact for large-d. We stress, however, that a we are not 
concerned here with the calculation of 0(1/A^) corrections to the previous results. 
Our purpose is to extend the leading order calculations to off-critical systems. 
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2.2 Vector fields 



For non-conserved vector fields, the starting point is the TDGL equation in the 
form (^. This time a vector field m(r,t) is introduced, whose zeros define the 
locations of the topological defect cores. The transformation 0(m) is now defined 
by an equihbrium solution of , which satisfies the vector analogue of |^, |^ 

Wl$=dV/d$, (28) 

where = YJ^=id'^ /d'^nlt is the Laplacian in fh space. We look for a radially 
symmetric solution of (^Sf ): 

— * 

cj){m) = g{p)m , (29) 
where p = |m| and m = m/|m|, with boundary conditions 

g{oo) = l , g{0) = 0, (30) 

(though any solution differing from ( P^ by a global rotation is equally acceptable). 
Thus, g{p) is the equilibrium profile function for a topological defect in the field 0, 
with \m\ representing the distance from the defect core |]^, |^. A solution of this 
form makes sense, of course, only for n < d when singular topological defects exist. 
Rewriting the TDGL equation @) for each component of in terms of m, and using 
(p8|) to eliminate dV/d(j)a, gives 

^ d(j)a drrib ^ d(j)a „2 , ^Va V7 V7 V72 , /oi N 

? d^,^ = ? + ? ^^^^^"^^-^"^^ - ^-'^^ ■ 

As in the scalar case, in order to establish a linear equation for rh one replaces the 
quadratic factor Vmj ■ Vrrij by its average (over the initial conditions): 

Vmb-Vnic {{Vm)^)6bc ■ (32) 

Here we used the fact that ((Vmf,)^) is independent of b, from global isotropy, to 
write it as ((Vm)^) where m is any component of fh. As before, one can attach an 
additional 'colour' index a (= 1, A^) to the vector m, such that the replacement 
(^) corresponds to taking the limit A^ oo in the theory. In this case, however, it 
also corresponds to the hmit n ^ oo. With (p2D, (E^ simplifies to 



E = E I^V^™. - (1 - {(V,nr-)) Vl^. . (33) 



dnib dt ^ dmh 

Finally, one would like to eliminate the explicit dependence of this equation on 
(j){m). Once again, BH exploit the expected insensitivity of the scaling functions to 
the details of the potential, by choosing the function 0(m) to satisfy 

V'J = -(m-V„) = -J2mb^, (34) 
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a direct generalization of (^. Substituting the radially symmetric form (pQ]) in (|34D 
gives 

CTl — 1 \ Tl — 1 

the equation for the profile function g{p) with boundary conditions (|30D , which 
generalizes (0). For small-p the solution is linear 

g{p) ~ p (P ^ 0) , (36) 

while for large-p the profile saturates as 

g{p) ^ 1 - (n - l)/2p2 (p ^ oo) , (37) 

from which we can take = ('^ ~ l)/2 as a definition of the defect core size. The 
potential V{(f)) corresponding to this choice of the profile function can be deduced 
from (|28|) , though it seems unlikely that a closed- form expression can be derived. 



With the choice (|3^) , equation (|33D now becomes 
/ dm 



-V^m-a{t)mlVm<l)a = , (38) 



where a{t) is given by ([T^) . Let us define the vector 

d = dm/dt-V'^m-a{t)m . (39) 

In principle, equation (|3^) allows for solutions where fl and Vm0a are orthogonal 
vectors without Q being zero. However, inserting the radial form (^) in (|38D gives 

(n-wj = ^n^ + g'ip)nn = 0, (40) 

where Q± and fiy are the components of Q perpendicular and parallel to rh, i.e. 

Qn = (n-m)m. (41) 



The coefficients g' and g/p in (^0]) only vanish for p — >■ oo, so the orthogonal vectors 
Q± and must vanish separately everywhere. Therefore = is the only physical 
solution of (pH). Then equation (p^) reduces to the self-consistent linear form (|13|)- 



(H), holding separately for each component of m. Taking gaussian initial conditions 
for each component, with zero mean and correlator (|T3|), the usual procedure leads 
to the same asymptotic form (2^) for the correlator (ma(l)ma(2)) (a = 1, ...,n), and 
thus to the same form ( PU] ) for the normalized correlator 

'{m{iy) (m(2)2) 
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The final step is to evaluate the pair correlation function, C(l, 2) = ^0(1) -0(2)^. 
This proceeds exactly as in the BPT treatment [0: from (|25| ) p scales as v^, then 
according to p9|)-(l30D we can replace the profile g{p) by 1 and the function 0(m) 
by rh at late-times. Since from (|13D the components of m are gaussian fields at all 
times, the calculation of the average C(l,2) = (m(l)-m(2)) leads to the standard 
BPT scaling function P (see also the appendix): 



C(l,2) 



n7 
2^ 



B 



n + l 1 
2 ' 2 



n 2 



1 1 n + 2 
2' 2' 2 



;7(i,2)' 



(43) 



where -B(x, ?/) is the beta function and 2-^1 [c^, ^; c; 2;] is the Hypergeometric function. 

Remarks similar to those of the previous section, regarding the possibility of 
systematically improving the results by expanding in apply here. In addition, 
the present (leading order) results are expected to become exact both for large-d 
and for large-n, when the linearization of the m equation, leading to the gaussian 
property of m, is justified. 



3 OfF-Critical Quenches 

In the last section we discussed in detail the systematic approach for the 0{n) model 
representing a system undergoing a 'critical quench'. A remarkable feature of this 
approach is that it can be extended in a simple manner to treat the situation when 
an external field is present and/or the initial state contains a bias. 

Regarding the external field, the key point of the approach is to select a con- 
venient effective driving force which simplifies the equation for m(r, t) and which 
captures the essential physics of the problem. To achieve this, we focus on the rele- 
vant effect of the external field in the asymptotic dynamics of the ordering system, 
which is different for scalar and vector systems. The initial bias, on the other hand, 
only enters the equation for m through the initial condition. To leading order in 
the systematic approach, we obtain a linear equation for the auxiliary field similar 
to the critical case, but with an extra driving force h for n = 1 and ^ {w?) h for 
n > 1. To evaluate the expectation values associated with the order parameter (in- 
cluding (^(Pj, which was zero in the critical case), we extend the gaussian calculations 
for a distribution centered about a non-zero average. In particular, we derive the 
off-critical extension of the BPT scaling function (|43|) . 

3.1 Scalar fields 

Our starting point for a non-conserved, scalar system in the presence of an external 
field, is the usual TDGL equation (|^), but now the potential V{(f)) is asymmetric. 
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The energy density due to an external field (e.g. a magnetic field) can be modelled 
by a linear term coupling the field to the order parameter, which biases the order 
parameter in the direction of the field. Hence the potential V{(f)) has the form 

V{4>-h) = Vo{^)-h{t)^ , (44) 

where Vo(0) is the usual symmetric, double- well potential, and h{t) is a general time- 
dependent external field. For a spin system, for example, the above linear form only 
holds if is far from saturation and h is not too large, otherwise the system will 
respond non-linearly to the field and terms such as 0{h(f)^) become important. In 
the present continuous order parameter system, the saturation values of (p simply get 
modified by h. For positive h and given t, the potential V{(f); h) has an asymmetric 
double- well form, with the right- (left-) hand minimum lower (higher) than in Vo(0). 
For our purposes, it is convenient to rewrite the potential in the generic form 

l^(0;/i) = K)(0)-M^)V^i(0) , (45) 

where Vi(0) is a monotonic, odd function of (p. Then the evolution equation for the 
order parameter reads 

^ = V'<P-V^i<P) + h{t)V;i<P) . (46) 

As in the treatment of the case h = 0, we can exploit the insensitivity of the 
domain growth to specific details of the potential by choosing an especially con- 
venient form for V{(f)). This relies on the physical idea that the interface motion 
only depends on the local curvature K and on the local field h. To motivate this 
idea, we derive the generalized version of Allen-Cahn's equation for the interface 
motion, valid for an arbitrary potential. If g is a coordinate normal to the interface 
(increasing in the direction of increasing 0), the TDGL equation (|^) can be recast 
near an interface as 

where the curvature K is the divergence of the 'outward' normal (i.e. pointing in 
the direction of increasing 0) to the interface. Multiplying through by {d(f)/dg)t 
and integrating over g through the interface gives the local velocity of the interface, 
V = {dg/dt)^ as 

v = -K + AV/S , (48) 

where AV is the change in the potential V{(j),h) across the interface, and S = 
/ dg{d(f)/dg)f is the surface tension. The essential point here, is that the interface 
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motion depends on the external field only through AV. So, as far as domain growth 
is concerned, all that matters is the difference between the minima of the potential 
V{(f); h) for the two bulk phases. This gives us a great deal of flexibility in the choice 
of Vi{(f)). In particular, we may choose the minima of V{(f); h) to remain at = ±1, 
avoiding inessential contributions to the magnetization from 'stretching' the fleld. 

As usual, we introduce the auxiliary variable m(r, t) whose zeros determine the 
wall positions. Rewriting equation (^) in terms of m yields 

0'm = + 0"(Vm)2 - ¥^{(1)) + h{t)V;{(j)) . (49) 

We deflne the transformation = 0(m) by the equilibrium wall-proflle for the 
asymmetric potential. This situation corresponds to a planar moving wall (moving, 
say, in the positive x-direction) with position xo{t) and velocity v{t) = xo(t). So, 
the equilibrium solution is of the form 

0(r,t) = 0(x-xo(t)) . (50) 

Substituting this in (^Bf ) gives the equation for the steady-state proflle 

V'{(f);h) = (j)"{m) + v{t)(j)'{m) (51) 

which generalizes (^, and has the same boundary conditions (|^). This immediately 
suggests the convenient choices 

VM = (52) 
K(0) = 0'(m) , (53) 

along with the identiflcation 

v{t) = -h{t) . (54) 

Note that (0) determines the proflle for a given Vo(0), or vice- versa, while ( |53D 
flxes Vi(0) for a given proflle 0(m). In addition, we adopt the BH choice (^) for the 
function (f){m), i.e. 0" = —mip'. Inserting (p^, ( |53D and (|^ in (^Dj), we obtain the 
simpler equation 

(l - (Vm)2) m + h{t) . (55) 

The condition (|^) gives once again the error function proflle (0), and amounts, 
via (|5^), to choosing the previous form (|10D for the symmetric potential Vo(0). The 
potential Vi(0) is obtained by integrating (|53|), with boundary condition Vi(0) = 0. 
This gives 

Vi{(f)) = dx {(f)'{x)) = — dxexp[—x'^) 
Jo n Jo ^ ^ 
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— =erf(m) (56) 



^erf fy2err^(0) 
a/tt ^ 



(57) 



where we used (H). Alternatively, from (^) and (H), we get 



\/i(0) = dx^2Voix) (58) 
Again, (|57|) or (|58|) only defines Vi(0) for 0^ < 1, but this is the only region we 



require for T = 0. In particular, we have Vi{(f)) ~ \j2/'K(f) for (jr <^ 1, while 
yi(0) ~ ± [l/v^- (1/4)(1 -02)y|in(l -02)|/2j for 0^ ^ i. Therefore, 1/i(0) 



has a sigmoid shape, with linear behaviour (^41) for small-0, while saturating to 
±l/y^ in the bulk and so acting only on the interfaces. We can relate our field 
h{t) to an effective field driving the interface dynamics by matching the differences 
= -2heff, which follows from (|4|), and AF = -h^Vi = -2h/y/TT. This yields: 
heff = h/y/iT. Note that if we had kept the form (|4^) throughout {V( = 1), the 



minima of V{(j); h) would be shifted by the field h, which is inconsistent with the 
use of the profile function (H), where 0^ < 1. This would also give a non-linear 
term h/(j)' ~ exp(m^/2) in (155|), amplifying the effect of the field in the bulk (and 
thus shifting the minima of V) and destroying the mathematical simplicity of the 
equations. 

Expectation values 

In order to solve equation (|55D, we use the replacement ([T^ ) for (Vm)^, which 
corresponds to taking the limit c? ^ oo or a large number of 'colours', A^. With this, 
(^) becomes the self-consistent equation for m{r,t): 

m = V'^m + a(t)m + h{t) , (59) 

with a{t) given by (|14D; a simple extension of (plsl). In addition, we allow for a 
uniform bias in the initial state, taking gaussian initial conditions for m with non- 
zero mean, and only short-ranged correlations 

(m(r,0)) = mo (60) 
(m(r,0)m(O,0))^ = (m(r, 0)m(0, 0)) - (m(r, 0))^ = A(5(r) . (61) 

The Fourier components of m(r, 0) still satisfy (mk(0)) = and (p!5|) for k 7^ 0. 
Solving equation (|59D in Fourier space, gives 

mk(t) = mk(0) exp (-kH + b{t)) + (27r)^5(fc) /* dt' h{t') , (62) 

^ ^ JO exp(o(r)) 
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with b{t) defined by ([T7|). Since ((Vm)^) is unaffected by the presence of a uniform 
bias and field, it is the same as for h = = itlq. Hence, the function h{t) is 
still determined by the self-consistency equation (0), leading once again to the 
asymptotic form (|T9|) for times t ^ Iq, with given by (^). Inserting this result 
in ( ]62| ) yields, for large-t, 

m^,it) = (t/to)'^ {mk(0) exp [-kh) + (27r)'^5(k) r]{t)} , (63) 

with 

i^it) = jjt' hit') (^) ' , (64) 

where the lower cut-off to was introduced into the integral to account for the break- 
down of the form used for h{t) at short-times. Although the form used for h{t) is 
only valid in the scaling regime {t ^ to) we shall be taking t > for computational 
convenience. This gives an inaccurate description of the initial transient behaviour 
prior to the scaling regime, but is irrelevant for the scaling behaviour (for more 
details see the discussion following (PT])). 

Averaging ( p^D over initial conditions and taking the Fourier transform, gives 
the average value of m(r, t) at late-times 

(m(t)) = {t/to)^[mo + r]{t)] . (65) 



From ( p3| ) we also obtain the previous expressions, ( pi] ) and (|25|), for the second 
cumulants of m, i.e. the connected pair correlator and one-point correlator 

Co(l,2) = (m(l)m(2))^ = (m(l)m(2)) - (m(l)) (m(2)) , (66) 
S'o(t) = (m^)^ = (m^) - {mf . (67) 

This immediately gives the same form (^6D for the off-critical normalized correlator 

^{in{ir}^ {m(2)=>, Wi+hf) \ 4(ii+t2)/ 

For equal-times (ti = t2 = t), (^) simplifies to 

7(1, 2)= exp (^-^^ . (69) 

As usual, at late-times we can replace the profile function 0(m) (eq. (|)) by 
sgn(m) to evaluate the expectation values associated with the field 0. According 
to (pU[), m(r, t) is a gaussian field at all times, with probability distribution defined 
by the non-zero mean (^) and the second cumulants (p^)-(p7|). Evaluating the 
gaussian expectation values we obtain scaling functions with arguments (|65|) , 
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and (|68D. The details of these calculations are given in the Appendix, so we shall 
state the results only. 



Magnetization 

For the average value of the order parameter 0, i.e. the magnetization, we obtain 

= (sgn(m)) = erf(M(t)) , (70) 
where M{t) is a dimensionless average of m: 

M{t) ^ -P= (71) 



d/8to (^-j [mo + v{t)] . (72) 

Some results for {(f){t)) are presented in Figures 1 and 2. The result using the profile 
(I) is (0) = erf [m/(1 + CV2 (m^) JI/2I ^ ^jth ( = ^2 the wall-width for this profile. 
Clearly this reduces to ([70| ) in the scaling regime. In the saturation limit (|M| 00) 
we expand the error function in (^0|) to give 

exp(-M2(t))' 



((/.(t)) ~ sgn(M(t)) 



|M(t)|»l, (73) 



|M(t)|v/? 

and so, the magnetization saturates exponentially fast. Conversely, in the 'initial- 
growth' regime when M is small (e.g. for small mg and h{t) and t not too large), 
(ffOp reduces to 

(0(t)) ^ (2/v/^)M(t) , |M(t)|«l, (74) 

showing that in this limit M plays the role of a magnetization. 

The time-dependence of (0) is determined by the spatial dimension d. On one 
hand there is the prefactor t'^/^ in M. On the other hand note that M has two 
contributions with different physical origins: one comes from the initial bias while 
the other is due to the external field. Which term dominates at late-times {t ^ to) 
depends on the value of d and on the explicit form of h{t). It is also interesting to 
note that if we set h(t) to vanish over a time ti > t^, then 

M(t) = {t/t,f' M{t,) , (75) 

i.e. the 'magnetization' induced at subsequent times grows in the same manner as 
the one induced by the initial bias. To illustrate the nature of the interplay between 
the magnetization growth induced by the initial bias and that driven by the external 
field, let us consider a time-independent field h. In this case (p3) reduces to 



ft d+2 

r]{t) = h dt' {to/t')— (76) 

Jto 
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2-d 1 2-d 



Ahto\^{t/to)— -Ij ~ h{t/to)— , d<2 (77) 
hto ln(t/to) , d = 2 

Ahtoll - {to/t)^\ ~ /i , d>2. 



where A = A/\d — 2\. The asymptotic forms on the right stand for late-times t ^ to- 
For d < 2 the integral ( |76D is dominated by times of order t, so the continued 
influence of the field dominates over the initial bias in ( |72|) at late-times. For d > 2 
the integral is dominated by short-times, so ri{t) scales like the initial bias mo in 



(^). This means that the prevalent effect of the field is the growth induced at early- 
times, a similar situation to that described by (|75|). For a general time-dependent 
field h{t) the crossover between these two regimes will be different, but clearly the 
factor in ( |64|) diminishes the effect of the field as time progresses. In 

particular, a sinusoidal field is less effective in reversing the magnetization during 
later cycles. This situation is illustrated in Figure 1. 

Using linear response and scaling arguments BK [|ll| predicted the scaling form 
of the magnetization in the 'initial growth' regime when (0) ^ 1. Our expression 
([f^)-([76|) for (0) agrees closely with their scaling form. The essential modification is 
that f^^^ = L'^/^ gets replaced by L^, where X = d — X is the exponent characterizing 
the decay of the autocorrelation function. Hence, the crossover between the two 
regimes no longer occurs at d = 2 but at the dimension where A = 1 |[Tl|. In the 
systematic approach, as in the OJK theory, the exponent A is given by 

XsA = d/2 . (78) 

The BK scaling prediction provides a means of measuring A without having to 
measure correlations between fields at different times. The result ([7(]|), although 
approximate, provides an explicit expression for the magnetization, describing the 
complete time-range from the early scaling regime to final saturation ((0) ±1)- 

To conclude this discussion, we address the situation when the initial bias uiq 
and the constant field h have opposite signs. According to ( [77D and ( [72D the mag- 
netization is reversed after a time trev given by 

{trev/tof'^ = 1 + \mo\/\htoA\ , for < 2 , (79) 
trev/to = ew{\mo\/\hto\) , for (i = 2 . 

For d > 2, however, a sufficiently weak field, \h\ < \mQ\/AtQ, can never overpower 
the initial bias and the system evolves towards the metastable state corresponding to 
(0) = sgn(mo). This result is for a system at zero temperature (there is no noise in 
the equation of motion). At temperatures T > thermal fluctuations will eventually 
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drive the system towards stable equilibrium (corresponding to (0) = sgn{h)) for any 
d. Does the above statement agree with what we should expect for T = 0? We recall 
the two effects of the magnetic field. Initially it creates a bias which grows in the 
same manner as the bulk magnetization induced by mo (eq. ([75D); the net result is 
the sum of these two. In addition it produces a wall driving force which adds to the 
wall curvature force (which exists for d > 1) |T^. The net force dictates which of the 



two effects dominates at late-times. Within our approach the curvature dominates 
for d > 2 (but not for d < 2), in which case the magnetization is determined by the 
net bias growth, i.e. the balance between mo and h. We would expect the curvature 
to dominate for any d > 1 when the field h is small enough, but we cannot say 
precisely, within this qualitative argument, how small h should be compared to mo. 
The results for 1 < d < 2 are probably an artifact of the approach, which is in fact 
a large-ci treatment. Therefore, the behaviour found for > 2 is probably the only 
physical behaviour when d > 1 and h is very small. 

Pair correlations 

For the pair correlation function of the field (p, we obtain (see the Appendix) 



C(l,2) = (sgn(m(l))sgn(m(2))) 
2 nCi.s) dy 



with 



TT Jo \/\ — 

M2(1) + M2(2) - 2yM{l)M{2) 



exp(-r(y)) + (0(1)) (0(2)) , (80) 



1-2/2 

where Mil) = M{ti) is given by (^) and 7(1, 2) is given by (p8|). In the critical limit 
(mo = /i = 0) we have M = (0) = 0, hence the usual OJK form C = {2/tt) sin^^ 7 is 
recovered. For widely-separated times (t2/^i 00) or for infinitely separated points 
(r ^ 00), we have 7^0, and so C{1, 2) (0(1)) (0(2)). 
At equal-times (M(l) = M(2) = M(t)), (|D|) reduces to 



f7W dy f-2M'^^ 

IT Jo vi - y \^ + y 



Cix,M) = 3 / -^^r^expl^^) + (0)^ (82) 



2 ^1 dy ( -2hP \ 
1 - - / /T g exp -—— , 83 



where x = r/t^^"^ is the scaling variable and 7(x) given by (pU]). It is worth studying 
the limiting behaviour of this expression. For 7 ^ 1 (x ^ 0) one obtains, to leading 
order in x, 

C{x,M) = l-exp(^-M^{t)) [{l/n)x + 0{x^)] , x^O. (84) 
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Note, that the term exp(— M^(t))x is, up to constants, Pdefit)r, where pdef is the 
wall density, given by (pSf ) below. This shows that the correlation function C {x, M) 
has a small- distance Porod's regime, as expected. In the limit 7 — > (x ^ 00), we 
take r — s> 00 keeping t fixed in order to have M{t) fixed. Then, Taylor expanding 
(p^) we obtain 

C(x,M) = ^ exp (-2M2(t)) {-f + M{tf-f^ + + {(f)f , x 00 . (85) 

In the saturation limit, when M(t)^ 00 and (0)^ ^ 1, C(a;,t) saturates exponen- 
tially to 1 in (^)-(|85|). The equal-time pair correlation scaling function C{x,M) is 



shown in Figure 3, for a constant external field. 
Density of walls 

We are also interested in the average density of wall 'area' per unit volume of the 
system, pdefif)- There are a number of ways of calculating pdef within the gaussian 
approximation. For example, pdef = ((V0)^^ /S = — (V^C)^^^/S, where S is 
the surface tension, can be evaluated from (|55D. The results differ by a constant 
amplitude, as a consequence of the inaccuracy of the approximation, but our main 
interest is to obtain the time- dependence of the wall density. A simple procedure 



17| , [18| which easily generalizes for vector fields, is to write the wall density at each 



point r in terms of the auxiliary field m, as 

Prfe/(r,t) = 5(m(r,t)) |Vm(r,t)| , (86) 

where |Vm| is the Jacobian between m and the spatial coordinate r. If we adopt 
the usual definition of m as a normal coordinate near defects, then |Vm| = 1 at 
defects, giving simply p(r, t) = 5(m(r, t)) Jl^ (note that this is true for the exact m, 
but not for the gaussian m, for which |Vm| fluctuates). Now we make the gaussian 
approximation for m, with the one-point distribution function 

Then, the average wall density is given by 

Pdefit) = (5(m)) = P(0) 

= — 7t — ' ^ ^ 

where we used (pSj), and M(t) is given by (|72]). This result agrees, up to constants, 
with an earlier prediction by OJK When rriQ = h = 0, ( PE| ) gives the usual scaling 
form p ~ L'^^^ I L'^ ~ The ordering process is much faster in the presence of an 
external field or initial bias, hence the exponential decay of the wall density. 
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3.2 Vector fields 



The evolution of a non-conserved vector system in the presence of an external field 
is given by the TDGL equation (j^). The usual form of the potential V{(j)) (for 'soft 
spins') is a direct generalization of (j^, i.e. 

V{$;h)=Vo{^)-h{t)-$, (89) 

where Vo(0) is the usual 'Mexican hat' potential, with a symmetric ground-state 
manifold 0^ = 1, and h{t) is a general time-dependent external field. The full 
potential V{(f)] h) has an asymmetric (tilted) 'Mexican hat' shape, with a single 
global minimum where (f) is parallel to h and has maximum length (among the all 
minima), and a metastable state where is antiparallel to h and has minimum 
length. These saturation lengths depend on the value of \h\. For the remaining 
directions of (p there is only a minimum along each direction, corresponding to a 
local saturation length, which varies continuously between the two above values. 

As usual for vector systems, we want to introduce a vector auxiliary field m 
whose zeros define the locations of the topological defect cores. By analogy with 
the previous cases, in order to define the transformation 0(m) let us imagine how 
a steady-state configuration of the field looks like in the presence of a slowly 
varying field h. For definiteness consider a system with n = d, which has point 
topological defects (i.e. vortices for d = 2 and monopoles for d = 3), along with anti- 
defects. These attract each other and eventually annihilate at late-times. Without 
an external field the defects are radially symmetric: at the defect core has zero 
length and undefined direction, while away from the defect core rapidly saturates 
in length. The effect of the field h is to deform the defects: all 'force lines' (except 
one) emerging from the defect core tend to line up with h away from the core. Close 
enough to the defect core, however, the field will still be isotropic, so we may still 
use the 'radially symmetric' form (^), i.e. 

0(rn) = g(p) m , (90) 

where p = \m\ and g{p) is a function yet to be determined. In the vicinity of a core, 
where m is isotropic, p and g{p) play the same role as in the case h = 0, while in 
the bulk the direction of rh is (on average) parallel to h. Note that, by symmetry, 
there is no direct force on the defect core due to the external field. 

Contrary to the scalar case, we conclude that an important effect of the external 
field occurs in the bulk, where the relevant phenomenon is the orientational ordering 
of along the direction of h. In the bulk we expect to rapidly relax to its saturation 
length, i.e. to the 'local' minimum of V{(j); h) which depends on the local angle 
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between cj) and h. In other words, since field 'stretching' is a local process while 
field rotation involves long-range correlations, the longitudinal effect of h (along (p) 
should be instantaneous as compared to transverse effects. The same assumption 



has motivated large-n studies of this problem using a 'hard-spin' model |]T9| , |TT 
This leads us to interpret g{p) in (^) as the stationary profile function of a defect 
in the field 0, with \m\ representing the distance from the defect along a (time- 
dependent) 'force line' defined by m. Hence, as far as bulk rotation and the slow 
defect dynamics are concerned, we may once again (but for reasons different from 
the scalar case) ignore the inessential 'stretching' of the field, by choosing the local 
minima of V(0; h) to remain at 0^ = 1. This amounts to taking the same boundary 
conditions (^) for g{p). 

With these insights, let us rewrite the TDGL equation @ as 



^ = VV--^ + C/(^;A). (91) 

where t/(0; h) is the driving force due to the field h (which would be simply h with 
the potential (P^)). As usual, we expect the scaling functions to be insensitive to 
the details of the potential, or equivalently, of the driving force. These should only 
affect phenomena on the scale of the defect core size (, but not on the scale of the 
bulk characteristic length L{t). Thus, we can make convenient choices for Vq and U 
in order to simplify the equation of motion for m. 

As in the case h = 0, we define the relation between the function 0(m) and the 
potential Vo(0) by the equation (^81). Rewriting (|9lD in terms of m gives (pT]) but 
with an extra term As usual, we use the replacement (|3^) which is justified in 
the large-n limit, where ((Vm)^) still stands for any component of m. Note that 
(Vm)^ is unaffected by a uniform external field or bias, so the system is isotropic 
as far as the average ((Vm)^) is concerned. We also choose the function 0(m) to 
satisfy (|3^). This leads once again, via (pOD, to the equation ( p5D for the profile 
function g{p), with small and large-p behaviour (p^)-(p7|). It also corresponds to the 
same choice for the symmetric potential Vo(0) as in the case h = 0. Putting it all 
together and using definition (P^, the TDGL equation (^Tj) then becomes 

(n-Vm) $=U{$;h) . (92) 

The key step now, is to rewrite the driving force in the anisotropic form: 

U{f,h) = U^{g)h^ + U\\{g)h\\ , (93) 

and to use ( pOD to write the left-hand side of (p^) as in (|40|) , i.e. as 

(n-Vm) $=^n^ + g'ip)nii , (94) 
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where the indices y and _L refer to the components of the vectors parallel and per- 
pendicular to (f) (eq. (|4l|) ). Matching (|93D and (0) gives immediately 



9KP) 

VL\Am) = -^hii Uii . (95) 
9{p) 

To obtain a simple, 'isotropic' equation of motion for m, we adopt the following 
choice for the driving force: 

U± = gip), 

Uii = pg'ip), (96) 



This corresponds, via (p3D, (|36D and (|37|) , to having 



U ^ ph (p 0) , 

U ~ + {I- (h\\ - h±/2) (p ^ 00) , (97) 



in total agreement with the our requirements from the above discussion: U vanishes 
at the defect cores, while being isotropic in their vicinity; in the bulk U reduces to 
its transverse component, forcing to rotate and to line up with h. Just as in the 
scalar case, the driving force derives from a potential: f/(0; h) = dVi{(f); h)/d(j). To 
get an equation for Vi one writes U in terms of and the components of 
parallel and perpendicular to h. However, according to ( P5| ) and (^) one cannot 



deduce Vi without solving equation (^) for g{p). Nonetheless, all that matters is 
the driving force which is physically correct. 

Expectation values 

We now solve the equation for m and obtain its expectation values. With the choice 
(0) and using definition (^) of f2, the equation of motion (|95|) becomes 



dm/dt = V'^m + a{t)m + ph{t) . (98) 

In the spirit of the systematic approach, in order to obtain a linear equation for rh 
we replace by its average in p: 



p=Vm2^^(m2), (99) 

which would be exact in the large-n limit. For convenience we set the external field 
in the (1, 1, ...) direction 

Mt) = ^(i,i,...) . (100) 

n 
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Clearly, the scaling functions will not depend on this arbitrary choice. As in the 
scalar case we also consider an initial state with a uniform bias. For simplicity we 
restrict the bias to be parallel or antiparallel to the external field: 

(m(r,0)) = mo(l,l,...) . (101) 

From (0), ( |100| ) and ( |101| ) it follows that the different components of m have the 
same expectation values at all times. Hence (m) = (m) (1, 1, ...) and (m^) = n (m^), 
where m is any component of m. Inserting ( pQ] ) and ( |100| ) in (0), we obtain the 
self-consistent linear equation for the components of m(r,t): 

dm/dt = V'^m + a{t) m + c{t) h{t) (102) 



a{t) = l-((Vm)2' 
cit) = = ^Soit) + {mf , (103) 



which is the vector field version of (|59|) . 

The calculation of the expectation values of m proceeds as in the previous sub- 
section, the only difference is that c{t), or equivalently M(t) or ri{t), has to be deter- 
mined self-consistently whenever h 0. As usual we adopt a gaussian distribution 
with mean mo and second cumulant ( |6lD for the initial condition of each component 
of fh. Fourier transforming (|102|) , and noting that b{t) (defined by b{t) = a{t)) is 
still given by ([T9| ) for t > to, yields the same solution (|63|) for mk(t), but where rjit) 
is now given by 

d+2 



r^{t)= jjt'h{t')c{t') (^) ' , (104) 

rather than (|6^) . This leads to the same expressions (|65| ) , (|2^) and (|25|) for the single- 
component averages (m), (m(l)m(2))^ and (m^)^. Hence the normalized correlator 
7(1,2) is again given by (|68|) , and the dimensionless average M{t) is given by (|72D 
but with rj{t) obeying ( |104|) . Using ( |103| ), ([TTl) and (|25|) we may rewrite the driving 



force as: c{t)h{t) = ^At/d^l + 2M'^{t) h{t). Inserting this in (|104D , and then ([lOiD 
in ([72|) , gives the exphcit equation for M(t): 



V2M{t) = i^-^j m'o + dt'h{t') i^-jj ^1 + 2M\t' 



(105) 



with itlq = ^J d/ 4:1710/ ^Jt^. For h this is a self-consistent equation which has to 
be solved numerically. It yields a much faster growth for M{t) than in the scalar 
case (eq. ([72D), where the integral in ( |105|) is just (0), i.e. c(t) = 1. Alternatively, 
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we could write equations for c{t) or ri{t), but M{t) is the relevant quantity for the 
scaling functions. Equation ( |105|) may be recast in the differential form 



d_ f V2M \ 
Jt [l^J 



h{t) 



1 



( V2M V 



(106) 



with initial condition at t = tn 



t = tn 



(107) 



This is equivalent to the equation obtained by BK [|IT| in the large-n treatment of 
the dynamics of 0. 

As suggested by (|106|) , we find that M is a scaling function of t and h(t). To show 
this it is convenient to separate the cases h = const and h = h{t). First we consider 
a time- independent external field h. Without loss of generality we set h > 0, while 
allowing the initial bias to be mo > or mo < 0. Defining the scaling function and 
variables 



N{y,z) = V2M{t)/y'/\ 
y = ht , 

z = m'J{ht,Y" = rr^^lyT 



equation (|106|) then reads 



dy V 1/^/2 ^ 



with initial condition at j/o = /ito 



A^(?/o, z) = z 



(108) 
(109) 
(110) 

(111) 
(112) 



where y^ is a dimensionless measure of the magnetic field. Note that N only depends 
on the initial bias through the initial condition. The solutions of ( |111| ) for large-?/, 
and for small-?/ and z (with logarithms for = 4), are 



N{y,z) 
N{y,z) 



B exp(?/) , ?/ > 1 



z-Ay'^ 



(4-d)/4 



(113) 
(114) 



where A = A/ (4— rf) and 5 is a constant. Using \/2M = N y'^/^ gives the correspond- 
ing behaviour of M{t). Note that for a negative bias starts negative, but since 
from ( p, 1 1| ) its growth is always positive, N becomes large and positive for large-y 
{B > 0). This can also be seen from ( p.05| ): if M starts large, there is a feedback 
through the time integral which boosts the effect of the h field. Therefore, contrary 
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to the scalar case (sec. |3.1| ) the external field will always reverse the magnetization, 
however large \z\ may be. The inversion time cannot be estimated, however, without 
solving the full equation ( |111| ). For example, even if we set M{t) = in ( |105|) we 



cannot estimate the continuated effect of M{t') from to to t. But if the initial bias 
is very small, i.e. if \z\ is small, N will be small from the initial time to the inversion 
time. Then we can set = in ( |114| ), giving (cf. ([79|)) tinv/to = exp (|mo|/|/ito|) 
for (i = 4 and 

=y7^ + \AI^ , or {Unv/h)^ =l + \m'^\l\hto\A , for (115) 

The late-time exponential growth of M{t), which will lead to a very rapid sat- 
uration of (^<P), can be understood from the non-linear equation of motion (^): 
keeping the important terms at late-times and writing m = pin, yields m = h± 
and p = ph\y Hence p ~ exp(/;,||t). So, despite the large-n treatment of m we self- 
consistently obtain the correct asymptotics. The scaling regime (characterized by 
the scale L{t) ^ Q comprises both the 'initial growth' when M ^ 1 (eq. ( |114|) ), 
and the intermediate stage prior to saturation (when M ^ 1). So we want to solve 
( |1 1 1|) (numerically) to obtain the full solution for M{t). The above scaling form, 
not only simplifies the equations but also provides a more transparent picture of the 
combined role of the variables t, h and mo- Namely, 1/h and (htoY^'^y/to are the 
characteristic scales of time and bias. 

We now consider a time- dependent external field h{t). For h{t) > and h{t) = 
at isolated points only, we may define the following one-to-one mapping between the 
scaling variable and the time t: 

y(t) = f^dt'h{t') . (116) 







Let t{y) be the inverse function of y{t). By analogy with the previous case we 
introduce h^, a characteristic amplitude of h{t). Then, defining the scaling function 

N{y{t),z) = V2M{t)/{hotY/\ (117) 
z = m'J{hotoY'\ (118) 

and noting that dN/dt = h{t){dN/dy{t)), equation ([106D becomes 

ON 



— = ^(hot{y))-^/^ + m , (119) 



with initial condition ( |112| ) at y{t) = yo = hotQ. The solution of ( |119| ) for y and z 
both small is 

N{y,z) ^ z+ dy'ihM))-"'^ , 2/ « 1 (120) 
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while the large- ?/ solution is ( |113|) , but with y given by ( |116|) . When h{t) = const 



this formalism reduces to the previous one. The procedure to relate the scaling 
variables y2{t) and yi{t), concerning two different fields h2{t) and hi{t), is to replace 
h2{t') by h2{t{yi{t'))) in the time integral defining y2{t)- In general the function 
y{t) will not be invertible, however, so the scaling form is of limited practical use 
in this case. For example, the integral of h{t) = /iq exp(— cut) is easily inverted, but 
the integral of h{t) = ho{l — sin(a;t)) has to be inverted numerically. In practice. 



of course, one can always solve equation ( |106| ) for each given field h{t). From a 
theoretical point of view, however, the scaling variable y{t) is interesting: it shows 
that all that matters at a given time, is the sum of the field at all previous times. 

To evaluate the expectation values for the order parameter (p in the scaling 
regime, we make the usual replacement 0(m) rh. Given the gaussian initial con- 
ditions and the linearity of ( |102| ) (following from taking the large-n limit), the com- 
ponents of m(r, t) are independent gaussian variables with non-zero mean. Evaluat- 
ing the gaussian averages we obtain n-dependent scaling functions with arguments 
M{t), 7(1,2) and/or So(t), calculated above in the large-n limit. We recall that, as 
far as the determination of the m moments and of the scaling functions is concerned, 
the leading order of the systematic approach, i.e. the large-n treatment of the rh 
equation, is equivalent to assuming that m(r, t) is a gaussian field. 

Magnetization 

Evaluating the gaussian average of the field 0, we obtain (see the Appendix) 



M(t)4= /^t^s(l-s2)Vexpf-M2(t)s2^ ^ n21) 



where M{t) = M(t)(l, 1, ...) follows from ([T00|)-([10TD, and M{t) is the solution of 



equation ( |106|) . This result reduces to the error function form (ffOD when n = 1. In 
the early scaling regime (M^ ^ 1) we can expand the exponential in ( |121|) . Then, 
from ( |105| ), or using (|1 20|) and ( |117|) , we obtain 



5 (1/2, (n + l)/2) 



to J J to 



(122) 



in close agreement with the predictions from the large-n treatment of |]TT| . In the 
saturation limit (M^ ^ 1) the integral in (|121| ) is dominated by small values of s 
and we can expand the binomial. Using ( |113| ) and ( |117D for h 0, and ( p,05| ) for 
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h = 0, this gives 

= l-^ + 0{l/M') , |M(t)|»l (123) 

exp (~2j'dt'h{t')) 

2n {t/toY'^ ^ ' 

n — l 1 

^ ^-—0 ^TTTTW ' ^ = 0- (125) 

2n tjiq [t/toY/^ 

Hence, the magnetization saturates exponentially fast in the presence of an external 
field, and as a power-law when there is only an initial bias. Again, this agrees 



with the large-ra limit of |n, [T^. As expected, the magnetization induced by the 
external field always dominates over the initial bias at late-times. In other words, 
the system loses the memory of its initial condition. This differs from the scalar 
case, where both h and itlq contribute to the exponential behaviour (cf. (^)). The 
magnetization curves for the full time-range are shown in Figure 2 for different values 
of n and h. 

Although the transverse magnetization is zero by symmetry, the coefficient (^—1) 
in (|123|) suggests that the leading correction to saturation is due to the part of 



transverse to h. This is easy to understand: recall that quickly saturates in length 
in any direction, hence the time increase in the net component of (p along h must 
come from rotation. 

Pair correlations 

The gaussian calculation of the two-time pair correlation function of (p is relatively 
straightforward, though the algebra is rather long (see the Appendix). The simplest 
way we found to express the result is the following: 



C(l,2) = (m(l)-m(2)) 

= ^fdufdv ^7 + M(l).M(2)-^^7 
TT 7o Jo 2 P{u,v) 



((l-u2)(l-i;2))— [ a{u,v) 



(126) 



a{u,v) = M\l)u'^ + M\2)v^ -2-f{l,2)M{l)-M{2)uh^ 

(3{u,v) = l-72(l,2)uV. (127) 

M{i) is the same as in ( |121| ), and 7(1, 2) is given by (pB]). The equal-time correlation 
function C (x, M) is show in Figure 3 for different values of n and constant field h 
(the same values as in Figure 2). 

An especially convenient algorithm for numerical integration is obtained with 
the substitutions u = cos(a) and v = cos(6): this yields a factor (sin(a) sin(6))" in 
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the integrand making its singularity for 7 = 1 (equal-times) explicitly integrable. 
The scaling function (|126|) generalizes the usual BPT scaling function to an off- 
critical system, and indeed it reduces to (^) when M = 0. It can be checked 
that any choice of axis, different from (|100|) , produces the same result ( |126| ). The 
transverse and longitudinal parts of C(l,2) (relative to h) can be read off from 
the big bracket in the integrand: the transverse contribution arises from the term 
(n — 1)7/2, while the longitudinal one arises from (1/2 — a//3)7 + M(l)-M(2). When 
n = 1 (|126| ) should reduce to (^) (with (0) given by ([70D). To show this one must 
as a single integral using appropriate variable changes. Although we 



express ( p^ 

have not been able to give an analytical proof, we have verified the equivalence to 
great numerical accuracy. Another check on (|126|) is the large-n limit, which gives 
(m(l) -771(2)) /((m2(i)) (m2(2)))i/2^ as expected. 

When 7^0 with fixed M(l) and M(2), i.e. r ^ 00 with fixed times, the 
correlations decay as (compare ( ^5]) for n = 1) 



C(l,2) 



0(2)) + a(ti,t2)7 + O(7') 



:i28) 



where 0(^1,^2) is a known function of M(l) and M(2). It also interesting to look at 
the decay of the autocorrelation function A{ti,t2) = C(0,ti,t2) when ^2 ^ = ^05 
/i = and mo ^ 0. In this case ( |105| ) and (|6|) give: M(l) ~ m^, M(2) ~ M(l)/7 
and 7 ~ {ti/t2Y^^- A careful expansion of (|126| ) then yields 



1 + 



[n 



l)B ft 



,2 
m5 



t2 > h 



(129) 



where B = (2'^'^+^)/^ — 1)/2. As expected, the saturation limit of the autocorrelations 
is the initial bias in the field (p. The coefficient {n — 1) indicates a contribution from 
rotation, as in ([L23|). 



Density of defects 

Finally, we evaluate the average density of defect core volume per unit volume of the 
system, pdefif). The gaussian calculation follows the same steps as in the previous 
section for a scalar system. This gives 



p,^j{t) = {6{m)) = (P(0))" 

exp (-M^(t) 



in 



fn/2 



(130) 



which generalizes (]88|) . P{m) is the one-point distribution for each component of m, 
given by (pTj) . For mo = /i = ( |130| ) gives the usual scaling form pdef ~ L'^~"/L'^ ~ 
1/L". 
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3.3 Comparison with experiment: 2d XY model with bias 

For systems with d = n = 2 the defects are vortices and anti-vortices interacting 
strongly below the transition temperature (known as 'Kosterlitz-Thouless' transition 
temperature; see e.g. [0). The energy of a vortex pair is of order ln{L/(), with 
( the vortex core size and L the vortex pair separation. If scaling holds, then the 
defect density Pdefit) should scale as L^^(t) (eq. ( |13(JD ). The product of these two 
results then gives the energy density, e ~ ln(L/^)/L^. 

The planar XY model is of particular interest since it appears to be a special 
case: there are no conclusive predictions or definitive measurements regarding the 
asymptotic growth in this systems, but there are indications that the growth is 
anomalous. Yurke et al. have suggested a slower growth law L ~ (t/lnt)^/^ 
for non-conserved systems, and indeed it has proven difficult to reach the scaling 



regime through numerical simulations (see e.g. [jT2|, ^ and references therein). More 
recently. Bray and Rutenberg p2[ concluded that the energy dissipation due to the 
ordering process occurs significantly on all scales between ( and the inter-vortex 
spacing L{t), suggesting that there may not be a single characteristic length scale 
in the system, i.e. that the scaling hypothesis may not hold in this case. On the 
assumption that scaling does hold however, the (t/lnt)^/^ growth is recovered p5| . 
Numerical simulations ||23| also provide evidence that the scale L{t) required to 
collapse the data for the pair correlation function, and the typical intervortex spacing 
p(t)^^/^, are not simply proportional to each other. 

The systematic approach cannot account for these possible logarithmic effects 
since, as a consequence of treating the dynamics of the field m in the large-n limit, 
a t^/^ growth is predicted for all n and d. For the same reason, it also gives the 
incorrect value d/2 for the exponent A (eq. ([78|) ). It is still interesting, however, to 
make a qualitative comparison between the theory and experimental or simulation 
data, namely to confront the effect of the initial bias in the results. 

Pargellis, Green and Yurke (PGY) have devised an experimental system 
exhibiting planar XF-model behaviour. The system consists of a nematic liquid 
crystal material placed between two plates. A suitable choice of the temperatures 
at the two plates creates a nematic-isotropic interface (parallel to the plates) near 
the centre of the cell, where the nematic director has a fixed angle with the normal 
axis. Hence the projection of the director onto the interface provides the xy degree of 
freedom. In addition, a normal alignment is imposed at the colder plate, creating a 
splay of the director from this plate to the interface. With these boundary conditions 
the symmetry (p —>■ —(f) is effectively broken (i.e. there are no defects with charge 
1/2) and at the interface the defects look like vortices. After inducing a thermal 



27 



quench in this system, PGY measured the density of defects p(t) = pdejit), the 
autocorrelation function A{t,to) and the magnetization (^(j){t)'^. To confirm that 
the system is governed by the XY-mode\ dynamics, PGY performed numerical 
simulations of this model using as initial condition, not a random field configuration, 
but rather one obtained from experiment. Although they found good agreement 
between the results, these did not exhibit the expected scaling behaviour p{t) ~ t^^ 
and A(t,to) ~ t^('^^^)/'^ ~ ^-o.585 (^^y^qyb the latter follows from Mazenko's theory 
0, ^), indicating that some aspect of the initial condition was preventing the 
system from reaching the scaling regime. To study the effect of a non-uniform 
distribution for the angle 6 between the initial field 0(to) and the x axis, PGY 
performed simulations using the probability distribution: 

P{e) = ^ [1 + A cos(A;^^)] , (131) 
Ztc 

where A is the 'bias parameter', and k = 1 corresponds to a dipolar deviation from 
uniformity. Higher order distortions {k > 2) were shown to produce a late-time 
behaviour indistinguishable from the one with zero-bias {A = 0), and therefore were 
excluded. The log-log plots of the simulation data (with different values of A) for 
p(t), A(t, to) and I (0(t)) I/I (0(to)/ I are shown in Figures 4, 6, and 8. The corre- 



sponding log-log plots from the systematic approach (sec. |37^ ), obtained using the 
relation M(to) = A/\/2ti (see below), are shown in Figures 5, 7, and 9. Taking 
into account that the scales of time are arbitrary, there appears to be a good qual- 
itative agreement between the theory and the simulations, namely in the manner 
in which the curves change with the bias parameter A. PGY found that the best 
fit of the simulations to the experimental data is for A = 0.07. At late times, the 
simulation data and the theoretical curves for p{t) (Figures 4 and 5) drop below 
the line t^^ (the expected decay for A = 0), as a consequence of the initial bias. 
This is what we expect from (|130|) and (|105|) , which give: p{t) ~ exp(— tmo)/t. 
From the zero-bias autocorrelations (expected to decay as t^(^~'*')/2) PGY measured 
{d — A)/2 = 0.543 ± 0.009, while the corresponding results from Mazenko's theory 
and from the systematic approach are {d — A)/2 = 0.585 (solid line in Figure 6) and 
((i — A)/2 = 0.5 (solid line in Figure 7). For A 0, the autocorrelation curves do not 
decay indefinitely, but have a saturation limit which is the initial bias 0o = I (0(^o)) I- 
This is what we expect from ( |129| ) which gives: ln{A{t, to)) ~ ln(0o) + l/(mot)°'^^ for 
t ^ to- The solid lines in Figures 8 and 9, with slopes A = 0.83 and A = 1, describe 
the initial growth of the magnetization (cf. ( |122[ )). The fits to the simulation data 
{A = 0.07) and to the experimental data, however, give A = 0.518 ± 0.021. The 
magnetization curves with the larger values of A indicate that the system is close 
to saturation, and therefore outside the scaling regime. This behaviour is described 



28 



by ( |125| ). Since the experimental and simulation results do not respect the usual 
scaling forms in part because of the bias, they are inconclusive about the existence 
of logarithmic factors in L{t). 

In order to compare the theory with the PGY data, we calculated the probability 
distribution of 9 corresponding to a gaussian distribution of the initial auxiliary field 
m(r, Iq). Since P{0) in ( |131| ) is maximum for 6' = 0, 6' is the angle between m = m(to) 
and the most probable direction mo = (m(to)), i-e. /i = cos(6') = 0(to) ■ <^0(to)) = 
rh-rho. Noting that P{m) is the product of two distributions of the form (|87|), using 
polar coordinates, exploiting the symmetry of cos(6') and making the substitution 
s = cos (6'), we obtain 



dm P{m)6{fi — m-rho) 



°° dmm 



2nS() 



exp 



+ ml\ ds 



25n 



exp 



mmo s 
So 



6{fj, — s) . 



Performing the integrals and noting that P{0) = -P(/i)| sin(6')|/2, gives 



P{9) 



exp i-M^) 
2n 



1 + Mov^cos(^) exp cos2(2^)j (1 + erf (Mq cos(^))) 



(132) 



where Mq = \Mq\ = m^/ \/2Sq and m^ = |mo| is the initial bias. Expanding (|132|) 
when Mq is small yields, to leading order. 



P{e) = — 1 + Mov^cos(e) 
2tt l 



(133) 



Comparing with ( |131| ) gives M(to) = Md/v^ = A/\/27r, which we used to calculate 
the theoretical plots for given values of A. 



4 Conclusions 

We have extended the systematic approach (to leading order) to investigate the 
asymptotic ordering dynamics of non-conserved quenched systems with an external 
field and an initial bias. An important difference relative to critical quenches is 
that, as the phase symmetry is broken the system evolves more rapidly towards 
final equilibrium. In particular, as the topological defects disappear more rapidly 
the scaling regime, characterized by the length-scale L(t) ~ v^, has a more limited 
duration. 

A key ingredient of the approach is to identify the important effect of the external 
field in the ordering process, which is a wall driving force for scalar systems and 
bulk rotation for vector systems. This allowed us to make convenient choices for 
the potential driving forces {V({(f);h) and U{(j)\h)) and to introduce a meaningful 
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auxiliary field. Treating the auxiliary field dynamics in the large-ci (or large-n) 
limit, we were then able to calculate the magnetization (eqs. ([70D, (|121|) ), the pair 
correlation function (eqs. (0), (|126|) , which extend the standard OJK and BPT 
scaling functions to off-critical quenches) and the density of defects (eqs. (|88|), (|130|) ), 
for scalar and for vector systems. 

The magnetization and the pair correlations are shown in Figures 1 to 3, for 
particular choices of n, (i, mo and h(t). These results have been studied in various 
limiting cases: eqs. ®-(0) for (0), eqs. (|12g)-(|T2D for (0), eqs. (|1-(|5]) for 
C(l,2) with n = 1 and eqs. (|T28D-(p9D for C{1,2) with n > 1. The respective 
behaviours have been discussed and are well understood and in accord with our 
expectations. We find, in particular, that the saturation of the magnetization and 
the decay of the defect density are exponentially fast. The initial-growth of the 
mag netization (eqs. (|73D- ([T22D ) agrees well with a previous predictions by BK ||Tl| 



based both on the large-n solution for and on general scaling arguments. We 
also compared our results for a biased 2d Xy-model with numerical simulation data 
1^ (sec. |3.3| and Figures 4 to 9), and found good qualitative agreement. These 



simulation data, on the other hand, had been found to agree well (for a certain 
value of the initial bias) with data from experiments on nematic liquid crystals [|12 



Finally, we recall that the auxiliary field formulation relies on the existence of 
stable topological defects, which govern the asymptotic dynamics of the system. 
One legitimate question, therefore, is whether these defects will ever form in the 
presence of an external field, or a bias, before the system orders completely. Clearly, 
the defects will form, and an intermediate scaling regime will occur (where L{t) is 
the typical inter-defect distance) if the field and the initial bias are not too strong. 
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6 Appendix: Gaussian Expectation Values 

In this Appendix we indicate how to perform the calculation of the gaussian expec- 
tation values which were stated throughout this paper. The standard results for 
h = mo = can be obtained as particular cases of the off-critical results. 

For a scalar field, = sgn(m), m is a gaussian variable with non-zero mean. 
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correlator 7 = (m(l)m(2))^ /(S'o(l)S'o(2))-^/^ and second moment Sq = {m?)^ (eq. 
(|67|)). We recall that M = (m) /v^2So (eq. (|7lD ). For vector fields, (f = m/\m\, we 
conveniently set the external field h and the initial bias rn-o in the direction (1, 1, ...) 
(eqs. ( |100|) -( p!0TD ). Hence, the components of m, generally represented by m, are 
independent gaussian variables with the same moments as above. We recall that 



M = (m) /v2So = {m) (1, 1, ...). We will need the following integral representations 

(134) 



sgn(m) 
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dz 

— exp 
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' izm 


in J-00 
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V2vrSo 


Jo \fz 



zm"^ 



2So 



(135) 



6.1 Magnetization 

Using the one-point gaussian distribution (^7|), it easy to obtain the following results: 

/ m ( zm^W Ml ( zrr?\\ 

First we consider a scalar field 0. Using the representation ( |134| ) and the result 
( |136| ), we can write the average value of as 

1 c°° dz ( z^\ 

U>) = (sgn(m)) = — / — exp izM . (139) 

in J-00 z \ 4 y 

Differentiating with respect to M, completing the squares and performing the inte- 
gral, gives d {(j)) /dM = (2/y^) exp(— M^). Finally, integrating with respect to M, 
with boundary condition (0(M = 0)) = 0, gives 

(0) = ^ / dxexp{-x'^) = erf(M) , (140) 

VTT Jo 



which is the same as (ff^). 

Now we consider a vector field 0. Using the representation ( |135|) and the results 
( |137|) - ([138D , and exploiting the rotational symmetry in the m space, we can write 
the average value of as 

I m \ (1, 1, •••) dz I m ( zvr?\ \ / / zvr?^ ^ " ^ 



ml / a/tt Jo \fz \ \/2Sn \ ISn I \ \ 2Si 



M r°° dz f zM' 



2 



exp n 



VttJo yi(z + l)— V 
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Using the substitution s = ^ z/{z + 1), we then obtain 

2 /-i 



M 



TT JO 



1 - s^) ' exp (-s^M' 



n 



which is the same as ( |121| ), and clearly reduces to ( p.4CI| ) when n = 1. 



(141) 



6.2 Pair correlations 

Using the gaussian joint distribution for the variables m(l) and m(2) (which may 
be the same components of m(l) and m(2)), one can prove the following results, 
with the notation mi = m(l), Si = ^0(1), Mi = M(l), and 5=1 — 7^, 
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(143) 
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with the definitions 



A12 
B12 



Mlzi + Mlz2 + 1(7)^1^2 
(zi + l)(2;2 + l)-7' , 



where r(7) is given by (^Tj). 

First we consider a scalar field 0. Using the representation (134) and the result 
( [L42|) , we can write the correlation function of (p as 

C(l,2) = (sgn(mi)sgn(m2)) 

1 dzidz2 (. zl + 4 + 2ziZ2j\ 

= J / exp tZiMi + IZ2M2 . 

vr J— 00 ^1^2 V / 

Differentiating with respect to 7, completing the squares and performing the in- 
tegrals, gives (9C(1, 2)757 = (2/7r) exp (— r(7)) / ^1 — 7^. Finally, integrating with 
respect to 7, with boundary condition C(7 = 0) = (</>(l)) (0(2)), gives 

dy 



C(l,2) 



TT JO ^1 — y 



^exp(-r(y)) + (0(1))(0(2)) 



which is the same as (PU]). For critical quenches, F = (0) = 0, and 
the OJK scaling function (2/7r) sin~^ 7. 



(145) 
reduces to 
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Now we consider a vector field Using the representation ( |135|) (with the sub- 
stitution z — > — 7^)) and the results (|143|) -( |T4^ ), and exploiting the rotational 
symmetry in the fh space, we can write the correlation function of (p as 
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Using the substitutions u = Zi/{zi + 1 — 7^) and u = Z2/{z2 + 1 — 7^), then gives 
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which is the same as ( |126| ), with a\2 and /?i2 given by ( |127D . For critical quenches, 
Ml = M2 = 0, and ([146D reduces to the BPT function (H). 
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Figure Captions 

Figure 1. Magnetization for d = 3 and n = 1, with a sinusoidal external field h{t) = 
0.01 sin(0.15(t/to — 1)) and an initial bias from mo = 0.003 to rriQ = 0.012. Note 
that (0(to)) > 0. 

Figure 2. Magnetization for d = 3 and n = 1,2, 3, with a constant external field 
h = 0.01, 0.02, 0.03 (bottom to top) and initial bias rriQ = 0. 

Figure 3. Equal-time pair correlation function for d = 3, n = 1,2,3 and at times 
t/to = 1,51, with a constant external field h = 0.01 and initial bias mo = 0. The 
behaviour at large-r is described by ( |128| ). 

Figure 4. Density of defects for the 2d XY model. Simulation data from Pargellis 
et al. There is no applied field and A is the bias parameter (eq. ( |131|) ). The solid 
line with slope —1 is the expected asymptotic behaviour. 

Figure 5. Density of defects for the 2d XY model, as predicted by the systematic ap- 
proach. There is no applied field, and the initial bias is given by M(to) = A/(27r)^/^. 

Figure 6. Autocorrelation function for the 2d XY model. Simulation data from 
Pargellis et al |]12|. The solid lines with slopes —0.5 and —0.586 are the asymptotic 
decays predicted by the systematic approach and by Mazenko's approach ||. 

Figure 7. Autocorrelation function for the 2d XY model, as predicted by the sys- 
tematic approach. 

Figure 8. Magnetization for the 2d XY model. Simulation data from Pargellis et al. 



Iir^ The solid lines with slopes 1 and 0.83 describe the initial growth as predicted 
by the systematic approach and by Mazenko's approach |^. 

Figure 9. Magnetization for the 2d XY model, as predicted by the systematic ap- 
proach. 
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